Correlations in orbital angular momentum of spatially entangled paired photons 

generated in parametric downconversion 
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What are the orbital angular momentum correlations between spatially entangled photon pairs 
generated in spontaneous parametric down-conversion? We show that the answer to this question 
can be given in two alternative, although complementary, ways. The answer posed in this letter 
explains satisfactorily the seemingly contradictory results obtained in different experiments, and 
theoretical approaches. 
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Within the paraxial quantum optics regime, the orbital 
angular momentum (0AM) provides an useful descrip- 
tion of the spatial degree of freedom of photons. Photons 
whose spatial waveform contains an azimuthal phase de- 
pendence of the form ~ exp(irrnp), carry an 0AM of 
mh per photon Photons with diverse spatial wave- 
forms can be easily generated, detected, and controlled. 
Therefore, the 0AM offers a physical resource where to 
explore deeper quantum features not present in the 2- 
dimensional Hilbert space addressed when using the po- 
larization 0, Q. Indeed, it allows to readily tailor the 
number of effective dimensions of the Hilbert space Q . 

During the last few years, several quantum features 
based on the capacity of the 0AM of photons to go be- 
yond a 2-D Hilbert space have been demonstrated (see 
[5| and references inside) using spontaneous parametric 
down-conversion (SPDC). These include the demonstra- 
tion of the violation of bipartite, three dimensional Bell 
inequalities Q , the implementation of the so called quan- 
tum coin tossing protocol with qutrits , and the gener- 
ation of quantum states in ultra-high dimensional spaces 

i- 

All of these experiments make use of the existence of 
specific quantum 0AM correlations between the two en- 
tangled photonsgenerated in the SPDC process. Several 
experiments 0, UllOj [HI seems to support the validity of 
the selection rule rup = mi +m.2, where mph is the 0AM 
per photon of the classical pump beam, and mi and m2 
are the winding numbers of the modes into which the 
quantum state of the signal and idler photons are pro- 
jected, respectively. Some other experiments [T2|. [l3l fl^ . 
while not directly measuring the 0AM of the downcon- 
verted photons, demonstrate the existence of ellipticity 
of the spatial waveform , which should make possible the 
detection of photons with ^ mi + m2- Under some 
restrictive conditions, the selection rule m^ = mi -\- m2 
can be derived from first principles 0, [H, although, 
as it will be shown below, the same rule address differ- 



ent physical quantities. The presence of Foynting vector 
walk-off can also strongly modify 0AM correlations [l3l ■ 

All this raises the question of what are the 0AM corre- 
lations between the downconverted photons generated in 
SPDC, i.e., under which conditions the 0AM of the en- 
tangled photons fulfill the selection rule = mi -I- m2. 
Here we show that this question can be formulated in 
two complementary scenarios, so that in each scenario 
the sough-after 0AM correlations can be different. The 
existence of previous apparently contradictory results is 
due to the fact that the sought-after 0AM correlations 
are different. 

In one scenario, the spatial properties of all the pairs 
of photons generated are considered, therefore the global 
mode function is obtained adding coherently all such pos- 
sibilities. In another scenario, which is relevant for cur- 
rent experimental applications, a small section of the full 
down-conversion cone is considered. Only certain proba- 
bility amplitudes are now considered. Under these condi- 
tions, the non-coUinear SPDC geometry and the presence 
of spatial walk-off can greatly modify the 0AM correla- 
tions observed. 

We consider a nonlinear crystal of length L, illumi- 
nated by a monochromatic laser pump beam propagat- 
ing in the z direction, with frequency LUp. The spatial 
shape of the pump beam at the input face of the nonlin- 
ear crystal (z = 0), in the transverse wavevector domain, 
writes E+ (p) = Eq {px + ipyT" exp (-Ippwg/d), which 
corresponds to a beam which carries an 0AM of mph per 
photon. Eq is a normalizing constant, p (px,Py) is the 
transverse wavevector and wq is the beam width. The 
signal and idler photons are assumed to be monochro- 



matic, with 



(jjp/2, where LUs^i are the frequen- 



* Electronic address; |juan.perez@icfo.es| 



cies of the signal and idler photons. This is justified by 
the use of narrow-band interference filters in front of the 
detectors. 

The photons are known to be generated into cones, 
whose shape is determined by the phase matching con- 
ditions inside the crystal. For the sake of clarity, let us 
consider first noncritical, i.e. negligible walk-off, non- 
coUinear SPDC in a periodically poled nonlinear crystal. 
The angle of the down conversion cone is assumed to be 



2 



small, so that the polarization [l8| and refractive index do 
not show noticeable changes with the direction of propa- 
gation. Similarly, the nonlinear coefficient is assumed to 
be constant along the down-conversion cone. 

The SPDC process can be described in the interac- 
tion picture by an effective Hamiltonian given by [l9l | 
Hi = eo JydVx'-^^E+E-E- + c.c, where E' {x,z,t) oc 
J dKsdP exp {—iP ■ x — iKgZ + iujst) {Ks,P) refers to 
the negative-frequency part of the signal electric field op- 
erators. Similarly for the idler photon. The two-photon 
quantum state j^*) at the output face of the nonlinear 
crystal, within the first order perturbation theory, can 
be written as |*) = / dPdQ $ (P, Q) al (P) a| (Q) |0, 0), 
where P and Q arc the transverse wavevector for the sig- 
nal and the idler respectively, and the mode function $ 
is given by 



(a) 



$ (P, Q) ^Ep{P + Q) sine 



exp < I 



.A,L 



I [Ks (P) + K, (Q)] L 



(1) 



where Aj. ^ Kj,{P 



Q) - Ks (P) - A'i (Q), the wavevec- 
\u,n,/cf-\pf"^ 



tors {j = s,i,p) write Kj (P) 

depends on the modulus of the corresponding transverse 
wavevectors, and rij are the corresponding refractive in- 
dex. The mode function of the biphoton in the spatial 
domain (xi,X2) is the spatial Fourier transform of the 
mode function given by Eq. ([T]). 

We can write |P 4- Qp = + + 2psPi cos {(ps - ifi), 
where ps ~ |P|, and ips = ta,n~^ Py/P^ are the 
modulus and phase of the transverse wave vector P 
in cylindrical coordinates. For the idler photon we 
have, similarly, pi — |Q| and ipi = ta,n~^ Qy/Q^. 
One can write sine exp (iAfcL/2) = 

J2hL-oo T^i (Ps: Pi) exp {il {(ps - 'Pi)}- Ks,i depends 
on the moduli ps^i, respectively. The pump beam can 
also be written as 



E.p = Eq exp < - 



[pI + p1 + PsPi cos {ip, - (f,)] 



^ ^ 1^ TTip j ^,^mp exp{ilips + iimp-l)ip,} (2) 

Therefore, the mode function given by Eq. ([1]) can be 
written as 

oo 

$(P,Q)= ^ gm{ps,pi)exp[irrups + i{mp~m)ipi] 

rn—~OQ 

(3) 

The main conclusion to be drawn from Eq. ([3]) is that, 
if polarization, refractive index and nonlinear coefhcient 
show negligible azimuthal variations around the down- 
conversion cone, the 0AM correlations of the spatial 
waveform of the biphoton state fulfill nip = mi + m2 [l5| . 
Importantly, this result requires considering the whole 




FIG. 1: Schematic diagram of a non-collinear SPDC. (a) 
Top view of the non-collinear configuration, (b) The down- 
conversion cone. 



spatial waveform of the downconverted photons, i. e. the 
full down-conversion cone. Notwithstanding, these are 
not the 0AM correlations that typical quantum informa- 
tion experiments based on spatial entanglement measure. 
An experiment aimed at detecting the global 0AM of 
the downconverted photons is a significant experimental 
challenge that it is yet to be solved. 

All relevant experiments reported to date detect only 
a small section of the full down conversion cone. In 
other words, the wave vectors of the signal and idler 
photons belong to a narrow bundle around the corre- 
sponding central wavevectors, i.e. P = Pq -I- AP and 
Q = Qo + AQ. As shown in Fig. 1, the signal pho- 
ton propagate along the direction zi with longitudinal 
wavevector fc, (p) - o _ i„|21 1/2 

verse wavevector p 



[(w^ns/c) 2 — IpI^J and trans- 
= {px,Py), so that APj = and 
APy = cosOiPy — sin Oiks- And similarly for the idler 
photon, which propagates in the direction Z2 with longi- 
tudinal wavevector ki (q), and transverse wavevector q, 
so that AQx — qx and AQy = cos62qy — sm92ki. We 
restrict ourselves to the case Oi ~ —62 ~ . 

The quantum state of the biphoton at zi ~ L/ cos 9 can 
be written as = J dpdq$ ( p,q ) a+ (p) a+ (q) |0, 0), 
where the mode function writes 12011 



$ (p, q) = Ep {p^ H- qx,5Q) sine f 

( S L L 

exp \ i4r -I- i [ks (p) + ki (q)] 

2 cos d 



(4) 



where 5k = kp — [kg + ki) cosO — {py — qy) smO and 5q = 
{Py + qy) cos 6 - {kg + ki) sin 0. 

The mode function given by Eq. ^ shows ellipticity in 
the (p, n) domain, as has been demonstrated experimen- 
tally [13, [3 • An increasing degree of ellipticity of the 
spatial mode function enhances the quantum probability 
amplitude of paired photons with nip ^ mi -|-TO2- To get 
further insight in the nature of the 0AM correlations, let 
us consider that the idler photon is projected into a gaus- 
sian mode (m2 = 0), so that the quantum state of the 
signal photon is described by the reduced mode function 
$s (p) (X / c?q<i> (p, q) exp (— |qpwJ/4), where wi is the 
beam width of the idler mode. To elucidate the 0AM 
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FIG. 2: Orbital angular momentum of the signal photon in 
non-coUinear SPDC in a L = 10 mm long PPKTP crystal. 
The pump beam is gaussian {nip = 0) . Both photons traverse 
2-f systems. The idler photon is detected with q = 0. (a) 
Total weight of the 0AM modes with mi 7^ 0, as a function 
of the pump beam width (wo). The label designate the non- 
coUinear angle, (b) and (c) OAM mode distribution for ip — 
1° for two values of the pump beam width: wq = lOO/im, and 
Wo = lOOO/xm. 



content of the signal photon, one has to project the spa- 
tial mode function into spiral harmonics exp (irrup) . The 
weights for each m of such decomposition gives us the 
sought-after OAM decomposition [3]. 

Fig. 2(a) shows the total weight corresponding to 
OAM modes with mi ^ 0, which is a measure of the de- 
gree of violation of the selection rule mp = mi + 7712 , for 
different values of the non-coUinear angle and the pump 
beam width. Notice that the larger the non-coUinear an- 
gle 0, and the smaller the pump beam width wq, the 
larger is the probability to detect signal photons with 
TOi ^ 0. Figs. 2(b) and 2(c) show two OAM distribu- 
tions for 9 ~ 1° , one which shows an OAM distribution 
with modes with mi ^ 0, while the other shows a single 
peak for mi = 0. 

The strength of the violation of the selection rule 
•nip = mi -|- mo can be quantified through the non- 
collinear length [131 Lnc = wq/ simp. If the crystal length 
is much smaller than the non-coUinear length (L <C inc), 
the ellipticity of the mode function is small, and thus the 
selection rule nip ~ mi + m2 is fulfilled. This turns out 
to be the case of nearly all of the experiments that make 
use of the OAM of photons @, 16., 1, 1, 1, ifl, [U . Typi- 
cal non-coUinear angles and crystal lengths used in these 
experiments are ~ 1 — 2° and L ~ 1 — 5 mm. For 
Wo = 500/im, one obtains Lnc — 15 — 30 mm, so that 
L < Lnc- On the contrary, if L > L„c, strong departures 
from the selection rule nip = mi -|-m2 are expected. This 



is the experimental configuration in [l^, [H, [13] , due to 
the use of a highly focused pump beam or longer crystals. 
For tuo = 90/ini and 9 = 2°, one has Lnc — 2.5 mm. 

According to Eq. ([3]), the OAM of the spatially en- 
tangled photons fulfill the relationship nip = mi -I- m2, 
while Fig. 2 shows that strong departures from this selec- 
tion rule can be observed if highly focused pump beams, 
highly non-coUinear configurations or longer nonlinear 
crystals are used. Actually, we are describing the same 
quantum process in two complementary scenarios. Fig. 
2 would give us the OAM correlations measured in typ- 
ical experiments that use the OAM as physical resource 
for quantum information, thus it is relevant for exper- 
imental configurations currently used. In this scenario, 
the fulfillment of the condition nip = mi + ni2 depends 
on the pump beam width and the non-coUinear angle, as 
dictated by the interplay between the non-coUinear and 
crystal lengths. 

On the other hand, Eq. ^ corresponds to a 
global view of the SPDC process, where the full down- 
conversion cone is considered. The question of angular 
momentum conservation balance in SPDC requires the 
simultaneous consideration of the angular momentum of 
the electronic spins and orbitals, the crystalline structure 
of the nonlinear crystal and of the electromagnetic field 
PH . The analysis presented here might be an important 
step towards clarifying how angular momentum is effec- 
tively conserved, since to evaluate conservation laws, one 
should take into account all probabilities amplitudes that 
contribute to the quantum process. 

Another important effect that might modify the OAM 
correlations is the presence of Poynting vector walk-off in 
some, or all, of the interacting waves. In type I and type 
II SPDC configurations, some of the interacting waves are 
extraordinary waves, thus show Poynting vector walk-off. 
For the sake of simplicity, let us consider a type I con- 
figuration, where only the pump beam presents spatial 
walk-off. An initially gaussian pump beam, at each po- 
sition z inside the nonlinear crystal can be decomposed 
into spiral harmonics as 

Ep {pp, z) = Eo exp (^!^ + I 

00 

X ^ Jn{zppta.npo)cxp{inLpp} (5) 

n — — 00 

where po is the spatial walk-off angle, ipp the azimuthal 
angle in cylindrical coordinates and J„ are Bessel func- 
tions of the first kind, and is the longitudinal wavevec- 
tor. As shown in Fig. 3, the OAM distribution of the 
pump beam increases its width with distance z. A gaus- 
sian beam incident from air into the nonlinear crystal is 
no longer a gaussian beam, i.e. the pump beam is no 
longer a wave with nip = 0. 

This effect becomes noticeable only if L > L.^, where 
the walk length writes L^ = wq/ tanpo- Fig. 3(a) shows 
the initial OAM distribution of a gaussian beam, which 
corresponds to a single peak with m^ = 0. Figs 3(b) and 
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FIG. 3; Orbital angular momentum distribution of the pump 
beam in different positions inside the nonlinear crystal, (a) 
2 = 0; (b) z — 5 mm and wo = lOO/im and (c) z — 5 mm and 
Wo = 1 mm. The walk-off angle is po = 5°. 

3(c) show the 0AM distribution at L ~ 5 mm for a walk- 
off angle of po = 5°. For highly focused pump beams, 
the walk off length can become smaller than the crystal 
length, thus the 0AM distribution no longer shows a 
single peak, as shown in Fig. 3(b). In this case, — 1 
mm. The walk off effect might be negligible for larger 
beam widths, as one can see in Fig. 3(c). Now, one 
obtains L^u ~ 12 mm. 

In the walking SPDC consider, the spatial waveform 
of the two-photon can be written as 

$ (P, Q) ^EpiP + Q) sine (AfcL/2) 

exp|i^+i[A-,(P)+/v,(Q)]L| (6) 

where Afc = fcp (P + Q) + (P^. + Q^) tanpo - h (P) - 
ki (Q). The phase matching function no longer depends 
on ips — ipi due to the presence of the spatial walk-off. 

The mode function of the two-photon state is a coher- 
ent sum of all the partial amplitudes due to the possi- 
ble generation of a pair of photons at position z inside 



the crystal [2^. Since the beam moves in the xz plane 
with angle po when propagating along z, the amplitude 
for each z shows increasing ellipticity, contributing dif- 
ferently to the total 0AM decomposition of the spatial 
waveform of the two-photon state. 

In conclusion, wc have shown that the elucidation of 
the 0AM correlations of entangled photons generated in 
SPDC can be addressed in two complementary scenar- 
ios, giving correspondingly different 0AM correlations 
between the photons. One scenario considers the quan- 
tum amplitudes in the whole down-conversion cone, while 
the other scenario, experimentally relevant for quantum 
information protocols that make use of the 0AM, pays 
attention to a small section of the down-conversion cone. 

The results obtained here are of great importance for 
other configurations where entangled paired photons are 
generated. This is the case of pairs of photons generated 
through two-photon Raman transitions in electromag- 
netic induced transparency schemes [23l |. Recently, the 
0AM correlations between stokes and anti-stokes pho- 
tons has been measured in trapped Rubidium cold atoms 
in a counter-propagating nearly coUincar geometry [2^ . 
The measured correlations comply with the selection rule 
nip — rric = mi — TO2, where TOc refers to the 0AM per 
photon of a laser control beam that counter propagates 
with the pump beam. The corresponding selection rule 
for counter propagating signal and idler photons [20| in 
SPDC would be rUp = mi — m2. We expect that a simi- 
lar experiment in a highly noncoUinear geometry, such a 
transverse emitting configuration [25j , would yield strong 
departures from the selection rule rUp — iric = mi — m2 
measured for a nearly coUinear geometry. 
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